Coulomb interaction at the metal-insulator critical point in graphene 
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We compute the renormalization group flow of the long-ranged electron-electron interaction at the 
Gross-Neveu quantum critical point between the semimetal and the excitonic insulator in graphene, 
perturbatively in the small parameter e = d— 1, with d as the spatial dimension. The O(e) correction 
to the usual beta-function makes the long-range interaction only more irrelevant at the critical than 
at the Gaussian fixed point. A weak long-range tail of the Coulomb interaction is found to be 
marginally irrelevant also in arbitrary dimension when the number of Dirac fermions is large. Its 
ultimate irrelevancy notwithstanding, it is shown that the metal-insulator transition may still be 
induced by increasing only the long-range tail of the Coulomb interaction. 



Increasing the strength of electron-electron interac- 
tions relative to the bandwidth is expected to trans- 
form graphene from its usual semi-metallic phase into the 
gapped Mott insulator. This quantum phase transition 
has been studied by a variety of numerical [l|, 0, H, 0] 
and analytical [1, 0, 0, @, El techniques, and it repre- 
sents a condensed matter analog of the particle physics 
phenomenon of chiral symmetry breaking. Graphene 
provides 2+1-dimensional Dirac fermions |1C| where the 
presence of two valleys and two spin states result in 
SU(A) chiral symmetry, which can be broken to SU(2) x 
SU(2) while still preserving parity and time reversal in- 
variance. Similar systems have been studied [ll||-03 as 
toy models of strong coupling behavior in quantum chro- 
modynamics and technicolor theories. In graphene, this 
results in breaking of the sublattice symmetry and gap- 
ping the electron spectrum. Generation and control of 
such a gap is of great importance to potential applica- 
tions in electronics flij . 

One issue of contention is the precise role in the mech- 
anism of the transition of the long-ranged ~ 1/r tail 
of Coulomb interaction, which remains unscreened when 
graphene is at the Dirac point. Whereas the initial ana- 
lytical calculations based on the Schwinger-Dyson equa- 
tions [1, HI found that the long-range interaction is cru- 
cial and leads to an essential singularity in the free en- 
ergy, recent numerical calculations [ij find only a regular 
second-order transition. The analytic expansion around 
the ex actly solvable three-dimensional limit of the the- 
ory [j| 0] also yields a regular critical point at which 
a weak long-range interaction may be shown to repre- 
sent a marginally irrelevant perturbation. In the present 
study we complement these results by showing that a 
weak unscreened ~ 1/r tail of the electron-electron re- 
pulsion remains a marginally irrelevant perturbation for 
an arbitrary number of Dirac fermion components N at 
the metal-insulator quantum critical point near one spa- 
tial dimension, as well as for large N in arbitrary spa- 
tial dimension. This is accomplished by computing the 
first correction to the beta-function of the long-range cou- 
pling constant, which we shall henceforth call the charge. 
This correction is 0(e) relative to the leading term, with 




FIG. 1: Schematic renormalization-group flow in the infrared 
for the metal-insulator transition in graphene. g is the (di- 
mensionless) short-range, and e 2 the long-range component 
of the Coulomb electron-electron interaction. GN and G are 
the critical Gross-Neveu and the Gausian fixed points, respec- 
tively. 

e = d — 1 , and its sign is such that the flow of the charge 
towards zero, although logarithmic in either case, be- 
comes faster at the quantum critical than at the Gaussian 
fixed point (Fig. 1). The renormalization-group flow also 
implies that the transition may be tuned by the charge 
alone, and can then occur even at a sub-critical value of 
the short-range interaction. 

We consider the simplest Gross-Neveu Lagrangian that 
should suffice to describe the quantum phase transition 
into the gapped excitonic (charge-density-wave) phase 
with increase of nearest-neighbor repulsion in graphene, 
tuned to be at the Dirac point 0, Q> 

L = <3> a 7o(<9o + ia )*a + v^ a jid^ a (1) 

-S(¥ a *a) 2 + alVl^V 

The short-range coupling g > is proportional to the 
nearest-neighbor repulsive interaction p}, and the four- 
component Dirac fermion ^ a is defined using the con- 
ventions of 0, Q. For generality we will assume an arbi- 
trary number of Dirac fermions N, so a = 1, 2, ...N. For 
physical spin-1/2 electrons N = 2. The Dirac matrices 
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satisfy the standard Clifford algebra in Euclidian space- 
time {7^,71/} = 28 hi,. For simplicity, we do not consider 
the transition to a spin-density wave phase 0, [§]• Our 
calculation can easily be modified to consider that case 
as well and the conclusion would be similar. 

Integration over the gauge field a introduces the in- 
stantaneous density-density long-range interaction 
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*a70*a(y,T), 



(2) 



in any spatial dimension d, provided that we define 

e ^ e 2 (47r) £ / 2 r( e /2), (3) 

and e = d—1. Note that as e — > the coupling e d defined 
this way diverges. This reflects the fact that the Fourier 
transform of ~ 1/r interaction in d = 1 is ~ ln(fc) (l5| . 
We will confine out attention to the region d > 1 in the 
following. 

Let us first outline our method and principal result. 
Since the inverse of the gauge-field propagator is a non- 
analytic function of the momenta for dimensions 1 < d < 
2, the coupling e d cannot renormalize and 
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where A is the ultraviolet cutoff. The instantaneous long- 
range interaction does not respect the Lorentz invariance, 
however, and consequently renormalizes the Fermi veloc- 
ity. At a momentum scale k <C A one expects the velocity 
to become 

v(k) =v+a d elHj) + b d elgMj)] n + 0(^,el9 2 ), (5) 

where a d and b d are d-dependent numerical coefficients, 
and g — gA e /ir is the dimensionless short-range coupling. 
The coefficient a d for d = 2 has been computed earlier 
0,0, Since v(k) is a physical observable it must 
be independent of the arbitrary cutoff A, i. e. 



dv{k) 
dA 



= 0, 



(6) 



at all momenta k. This condition of renormalizability of 
the field theory in Eq. (1) can be satisfied only if the 
power of the logarithm in the third term in Eq. (5) is 
n = 1, and if dg/dk — 0. Since the beta- function for the 
Gross-Neveu coupling, as it will be shown later, reads 

^ = eg+(l-^)ge 2 -2(2N-l)f + 0(g 2 e 2 ,g% (7) 

where e 2 = e 2 /irv is the dimensionless charge which mea- 
sures the strength of the long-range coupling relative to 
the Fermi velocity, the condition of renormalizability at 
e 2 = is satisfied only at the fixed points: 1) Gaussian, 
g G = 0, and 2) the critical g c = e/(2(2N - 1)) + 0(e 2 ). 



Eqs. (5) and (6) yield then the beta-function for the 
Fermi velocity 
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-e 2 d (a d + b d g*). 



(8) 



g* is one of the above two fixed points. We can then 
recast Eqs. (4) and (8) together as the equation for the 
dimensionless charge 
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= e\{^y/ 2 T{e/2)(a d + b d g* 
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Both of the coefficients a d and b d turn out to be O(e). 
Therefore the divergent ~ 1/e prefactor in the last equa- 
tion cancels. We thus find a regular expansion of the 
beta-function for the dimensionless charge in powers of e 
and e 2 
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+ 0(e z /N) +0(e b ). (10) 



There are two notable features of the last expression. 
First, the O(e) correction has the same sign as the lead- 
ing term. As a result, the long-range interaction is in fact 
more irrelevant in the infrared, although still marginally 
so, than at the Gaussian fixed point. Second, the correc- 
tion is 0(1/N). Therefore, for a large number of Dirac 
fermions, the Coulomb interaction is also marginally ir- 
relevant in all dimensions. Together with the previous 
analysis near three spatial dimensions 0| this means that 
a weak long-range tail of Coulomb interaction is an irrel- 
evant perturbation at the Gross-Neveu metal-insulator 
quantum critical point in all perturbatively accessible 
regimes of the theory. 

The beta-function for the short-range coupling in Eq. 
(7), besides the usual terms 0, H] obtains the contribu- 
tion from the long-range Coulomb interaction. The sign 
of this term is negative, so that there is a trajectory in 
the region g < g c that flows right into the critical point 
and separates the flows towards the semimetallic and the 
insulating ground states, as depicted in Fig. 1. The neg- 
ative sign is crucial to this result, which we may have 
expected on physical grounds: that an increase of the 
charge alone should take the semimetal towards the in- 
sulator. This situation may be analogous to the one in 
3 + 1-dimensional quantum electrodynamics, where the 
chiral-symmetry-breaking transition may also be tuned 
by the electromagnetic charge, but the critical behav- 
ior seems, on the other hand, to be controlled by the 
Nambu-Jona-Lasinio theory with only short-range inter- 
actions [H ED]. 

Let us now present the calculational details. The log- 
arithmic terms in Eq. (5) derive from the diagrams in 
Figs. (2)-(3). Diagrams in Fig. (3), however, do not 
contribute to the renormalization of the Fermi velocity. 
Since the fixed-point interaction is a contact interaction, 
the diagram in Fig. (3b) vanishes, whereas the diagram 
in Fig. (3c) is independent of the external momentum. 
The diagram in Fig. (3a), which is also the only dia- 
gram to this order which would be of order N, vanishes 
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FIG. 2: One- and two-loop diagrams that renormalize the 
Fermi velocity. The wavy line is the long-range interaction, 
and the dashed line is the contact fixed-point interaction. 



(b) 




as well because it is proportional to the factor TV70 = 
that accompanies the fermion loop. This implies that 
both the coefficients ad and bd as far as the number of 
Dirac fermions N is concerned are O(l). Since the fixed 
point interaction is g c ~ 1/N, it guarantees that the first 
correction to the Gaussian result is also 0(l/iV). 

One is therefore left with the diagrams in Fig. (2) to 
compute. Defining the self-energy from the Dirac fermion 
propagator G(k) as 



G (k) = ikojo + ivkiji + £(&), 



(11) 



where the d + 1-momentum k = (fco , k) , we can write 
both contributions together as 



S(fc) = ie 2 J ■ 



dq q a + k a 



r{-7o7a7o + (12) 



(2^+! W{q + kf 
2.9(7M7o7^7a7o + 707 (<?)}, 
where dq = dqod d q, and the integral 

j ( \ — f d P P^(P + l)u 

^ {q> J (2f)Wp 2 (p + (| ) 2 - 



(13) 



We have set the Fermi velocity in the last two equations 
to v = 1 for notational simplicity. Since the integral over 
the momentum q in the first term is only logarithmically 
divergent in the infrared, for the divergent part of the self- 
energy we need to include only the q — limit of 1^ (q) . 
In any d > 1 the integral I^(q) is infrared convergent 
and may be easily computed to be 



V(*) = ^f + o(<f 



(14) 



when eCl. On the other hand, 

<W7ju7o7f7a7o + 7o7q7^7o7m) = 2 (1 - d)7o7o7o, (15) 

where we analytically continued the identity 7^ = d to 
arbitrary spatial dimension d. The self-energy therefore 
simplifies into 



E(fc) = -ie^ 7 o7a7o(l + — ) 



dq q a + k a 



7T ' J (2TT) d + 1 \^(q + k) 2 ' 

(16) 

and becomes simply proportional to the leading contri- 
bution. The remaining integral may be performed in ar- 
bitrary dimension d > 1, and it yields 

(1 - d~ 1 )e : i , gA7 , ,A. , , s 

^ = rWw (1 + V )h( I^ (17) 
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FIG. 3: The remaining three two-loop diagrams that do not 
contribute to the Fermi velocity renormalization. 



In the limit d — > 1 the diverging part of the self-energy 
therefore becomes 



e 2 A 
£(fc) = —{1 + g)\n{-)ik ilt 
it k 



(18) 



Replacing the Gross-Neveu coupling g with its O(e) crit- 
ical value the last result may be recast as in Eq. (10). 

We turn now to the calculation of the beta-function 
of the short-range coupling in presence of the long-range 
Coulomb interaction. Since the Coulomb coupling vio- 
lates effective Lorentz invariance of the Gross-Neveu La- 
grangian, the engineering dimension of the short-range 
coupling constant receives an additional contribution 
from a non-trivial dynamical exponent z, and it is equal 
to e — (z — 1), The beta- function of the short-range cou- 
pling then becomes 
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(e-z+ l)g - 2{2N - l)g' + c d e z g + 0{g z e z , g J ), 

(19) 

where the third term arises from the diagrams shown in 
Fig. 4, and the term proportional to g 2 yields the usual 
beta-function for the Gross-Neveu coupling in the ab- 
sence of other interactions 0, In fact, only the di- 
agram in Fig. 4a renormalizes the short-range coupling, 
g — ► g + 8g, where 



5g = 



4tt 



e 2 9 



dq 



(2n) d + 1 {q + k) 2 \^ 



(20) 



This integral may be performed in arbitrary dimension 
d > 1, similarly to the one in Eq. (17), and in the limit 
d — > 1 its diverging part has the form 



5g 



■-e 5 m(- 



(21) 
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FIG. 4: Renormalization of the Gross-Neveu coupling due to 
the long-range Coulomb interaction. 



with C as the integration constant, T(a, b) = 
t a ~ 1 e~ t dt. Depending on the initial condition, Eq. 
([23| describes the infrared flow of the short-range cou- 
pling towards either zero or infinity. The trajectory end- 
ing exactly at the critical point as the charge flows to 
zero, g(e 2 — ► 0) = g c , separates two regions with differ- 
ent asymptotic behavior of the flow in the infrared, as 
in Fig. 1. The separatrix is obtained by setting C = 
in Eq. (|23[) . so that the solution near the Gross-Neveu 
critical point becomes 



yielding ca = —2/e in Eq. lfl9|l . The remaining dia- 
gram in Fig. 4b does not renormalize the Gross-Neveu 
coupling, but generates a new short-range coupling of 
the form ('5 707 'J) 2 , which is irrelevant close to the crit- 
ical point. An analogous situation arises in the (2 + 1)- 
dimensional quantum electrodynamics with the addi- 
tional Gross-Neveu interaction [23|. The irrelevance of 
the generated short-range interaction also agrees with 
the emergent Lorentz invariance in graphene [a]. Finally, 
identifying the right-hand side of Eq. (8) as being pre- 
cisely (z — l)e 2 [7| yields the beta-function of the Gross- 
Neveu coupling as in Eq. lff|. 

Remarkably, it is possible to obtain the explicit solu- 
tion of the flow-equations in Eqs. |(7|) and iflOj) . Since 
the long-range Coulomb interaction is marginally irrele- 
vant in the infrared, it is convenient to cast the two flow 
equations in the form 



"dg/de 2 =eg- 2(2N - l)g 2 - (2/e - l)ge 2 



(22) 



where we dropped the hats on the dimensionless cou- 
plings, and neglected the term ~ e in Eq. iflOj) . General 
solution of this differential equation reads 



ff(e 2 ) = g, 



(e/e 2 ) 



2\(2/e-l)„-e/e 



C + r(2/e,e/e 2 ) 



(23) 



5(e 2 



= 1 - 



-e 2 + 0(e 4 ). 



(24) 



The separatrix therefore lies in the region g < g c , and 
approaches the critical value of the Gross-Neveu coupling 
as the charge flows to zero linearly. 

To summarize, we have computed the beta-function 
governing the flow of the long-range tail of the Coulomb 
interaction at the Gross-Neveu metal-insulator quantum 
critical point in an expansion near one spatial dimension. 
The effect of the short-range Gross-Neveu interaction is 
to render the charge more (marginally) irrelevant than at 
the Gaussian fixed point. We discussed how the form of 
the flow diagram implies that the (irrelevant) charge is 
nevertheless a possible tuning parameter for the metal- 
insulator transition in graphene. 
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